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ABSTRACT 

We have made a Galactic extinction map of the Cygnus region with 5' spatial 
resolution. The selected area is 80° to 90° in the Galactic longitude and -4° to 8° 
in the Galactic latitude. The intensity at 140 /im is derived from the intensities 
at 60 and 100 /im of the IRAS data using the tight correlation between 60, 100, 
and 140 /xm found in the Galactic plane. The dust temperature and optical 
depth are calculated with 5' resolution from the 140 and 100 /xm intensity, and 
Ay is calculated from the optical depth. In the selected area, the mean dust 
temperature is 17 K , the minimum is 16 K, and the maximum is 30 K. The 
mean Ay is 6.5 mag, the minimum is 0.5 mag, and the maximum is 11 mag. The 
dust temperature distribution shows significant spatial variation on smaller scales 
down to 5'. Because the present study can trace the 5'-scale spatial variation of 
the extinction, it has an advantage over the previous studies, such as the one by 
Schlegel, Finkbeiner, & Davis, who used the COiJiiJ/DIRBE data to derive the 
dust temperature distribution with a spatial resolution of 1°. The difference of 
Ay between our map and Schlegel et al.'s is ± 3 mag. A new extinction map of 
the entire sky can be produced by applying the present method. 



Subject headings: dust, extinction — infrared: ISM — methods: data analysis 
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1. Introduction 

Radiation from cxtragalactic objects is generally reduced and attenuated the Galactic 
interstellar dust. As the dust extinction is high in the ultraviolet (UV) - optical wavelength 
region, the observed data at these wavelengths must be corrected for the Galactic extinc- 
tion along each hne of sight. For example, statistical analyses of galaxies using the Sloan 
Digital Sky Survey (e.g., York et al. 2000) require an accurate extinction map. Schlegel, 
Finkbeiner, & Davis (1998) (hereafter SFD98) have made an all-sky map of the Galactic 
extinction for this purpose by using the all-sky far-infrared maps of the Infrared Astronomy 
Satellite {IRAS) and the Diffuse Infrared Background Experiment (DIRBE) on the Cosmic 
Background Explorer [CO BE). 

To derive extinctions, the dust size is important; it The dust size is classified into 
three groups (e.g.. Desert et al. 1990): polycyclic aromatic hydrocarbon (PAH), very small 
grain (VSG), and large grain (LG). LG absorbs UV-optical photons and reradiates the 
absorbed energy in far infrared (FIR). The observed FIR emission is dominated by LGs (e.g., 
Sodroski et al. 1997). LG considered to be in equilibrium with the interstellar radiation field 
(ISRF) (e.g., Draine & Anderson 1985) and radiates thermal emission with an equilibrium 
temperature. Therefore, assuming a spectral emissivity index, the dust temperature and 
the optical thickness can be derived from the FIR spectral energy distribution (SED). We 
can derive the dust extinction in the UV and optical wavelength regions from the FIR SED 
if the dust optical depth of LG is proportional to the extinction of each grain component 
(PAH+VSG+LG). 

PAH and VSG, in contrast to LG, are heated stochastically by a single photon to high 
temperatures and cooled by radiation shorter than 60 /im. As a significant fraction of the 
60 fim intensity usually comes from VSGs, we cannot derive the LG temperature from the 
IRAS 60 and 100 fim band data (e.g. Desert et al. 1990). 

Furthermore, there are two all-sky FIR diffuse maps obtained by (IRAS) and DIRBE. 
IRAS has two FIR photometric bands of 60 and 100 jim, while DIRBE has four FIR photo- 
metric bands of 60, 100, 140, and 240 iim. The spatial resolution of DIRBE is 0°.7, whereas 
that of IRAS is 5'. 

SFD98 derived the dust temperature from DIRBE data at 100 and 240 //m, and pub- 
lished the all-sky extinction (reddening) map. They made an extinction map with 5' resolu- 
tion using the temperature derived from DIRBE and the 100 /im intensity of IRAS. Because 
the spatial resolution is limited by DIRBE, this extinction map may not trace small-scale 
variation as pointed out by Arce & Goodman (1999). 

Dobashi et al. (2005) made a visual extinction map with 6 ' resolution using the Digitized 
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Sky Survey (DSS). Since the star counting technique cannot measure the dust behind the 
stars, its accuracy can worsen in directions where the extinction is high, depending on many 
parameters such as the angular density of stars, their location from the molecular cloud, and 
the distribution of the distances of the dust. Therefore, this technique may be ineffective 
for the correction of the Galactic extinction towards extragalactic objects near the Galactic 
plane. 

To improve the maps, we assumed the color-color correlations between 60, 100, and 140 
Hm that were found by Hibi et al. (2006). They analyzed Zodi-Subtracted Mission Average 
(ZSMA) data of the Galactic plane within | 6 | < 5°, LMC, and SMC. They identified two 
groups having very tight correlations in the color-color diagram. These correlations are not 
produced with the model by Li & Draine (2002). Bot et al. (2009) also argued that the 
model by Draine & Li (2007) does not explain the Multiband Imaging Photometer for Spitzer 
(MIPS) and IRAS data. 

By using Hibi et al.'s correlation, we can obtain the dust temperature and the dust 
extinction with 5' resolution only from the IRAS data. To test the validity of this method 
of making a high-resolution dust extinction map, we have analyzed the Cygnus region. In 
this paper, we report the results of our analysis and compare them with those of SFD98 
and Dobashi et al. (2005). The data and the analysis method are shown in Section 2. The 
temperature and extinction maps are presented in Section 3. We discuss the results in 
Section 4. Section 5 presents to the summary. 

2. Data and Analysis 

In this section, the intensity at 140 iim is derived from the intensities at 60 and 100 /xm of 
the IRAS data, by using the tight correlation between 60, 100, and 140 /xm found by Hibi et al. 
(2006). To test this method, we selected a sample area in the Cygnus region. The selected 
area is 80° to 90° in the Galactic longitude and —4° to 8° in the Galactic latitude. The 
interplanetary dust (IPD) emission and the cosmic infrared background possibly contaminate 
the Galactic diffuse emission (e.g., Sodroski et al. 1997). However, as the Cygnus region is 
located at 60° in the echptic latitude, the subtraction of the IPD component has only a minor 
contribution to the uncertainty about our derivation of extinction discussed later. The cosmic 
infrared background emission is negligible, since the FIR intensity is high enough. 

In general, the interstellar space in the Galaxy is optically thin in the FIR wavelength 



-4 - 



region. The FIR intensity is written by 

4(A) = Tioo t,m X X B,{\ To). (1) 

In Equation (1), /,/(A) is the intensity at wavelength A, Tioo nm is the optical depth at 100 
/xm, Td is the dust temperature, and f3 is the spectral emissivity index. B,y{\,Td) is the 
Planck function. Hereafter, Iu{^) is described as /(A). /(A) is not color-corrected. The 
color-corrected intensity, /c(A), is described below. We adopted a one- zone model, i.e., the 
dust temperature is assumed to be constant along each hne of sight. We adopted jS — 2 
according to the previous study (Drain & Lee 1984, Hibi et al. 2006, SFD98) for comparison. 
Some observational results show that /3 changes depending on the physical status of the ISM 
(e.g.. Reach et al. 1995 and del Burgo et al. 2003). 

Because the published DIRBE all-bands intensities were derived by assuming a source 
spectrum, i^/(A) = constant, the color correction should be applied to the observed data 
with the photometric bands of DIRBE. The color correction factor Kx {P, T^) is determined 
from the spectral emissivity index and the dust temperature, if the spectrum is written 
by Equation (1). The conversion from the observed (uncorrected) intensity to the color- 
corrected intensity is given by Equation (2) with the color correction factor 

The color correction factor of DIRBE is quoted from the C05£'/DIRBE Explanatory Sup- 
plement (1995). The color correction was made for the intensity at 140 and 100 /xm, but not 
for that at 60 //m, because the SED shape around 60 //m cannot be described by Equation 
(1) (Hibi et al. 2006). The correction affects the intensities by a few percent and hence is 
not crucial for the degree of extinction derived later. 

The color-color diagram was made from the DIRBE data in the Cygnus region using 
the same procedure as Hibi et al. (2006). Hibi et al. (2006) indicated that there were two 
tight correlations in the Galactic plane. The major group is called the main-correlation, and 
the minor group the sub-correlation. The main-correlation and the sub-correlation are fitted 
with Equations (3) and (4), respectively. 



/c(140 /xm) / /(60 /xm) 



/c(100/xm) ■ V-^c(100/xm) 

/c(140 /xm) ( /(60 /im 



-0.78 



0-65 , , \ ■ (3) 



/c(100/xm) V^c(100/xm) 



-0.56 



0.93 / ^" ^ \ . (4) 
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Hibi et al. (2006) and Hirashita et al. (2007) indicated that the main-correlation represents 
the hnes of sight along which the ISRF is constant, and the sub-correlation represents the 
lines of sight along which the ISRF varies greatly. 

The color-color diagram of the Cygnus region is shown in Figure 1. In contrast to the 
Galactic plane, the data points satisfying 7(60 iJt,m)/Ic (100 jim) < 0.28 do not follow the 
main-correlation and other data points seem to follow the sub-correlation with a small offset. 
Therefore, we tried to fit the data points in two regions divided at 1(60 fim)/Ic (100 fim) = 
0.28. Least-squares fitting with two parameters of power and offset in the color-color diagram 
provided the following two relationships: 



Ic (140 /xm) ^ f I (60 ^m) \ _ 



-0.44 

(6) 



/c(100 jim) ' \Ic (100 urn) 



^c(140 /xm) ^ ^ / 7(60//m) \ 
7c(100 ^Jim) ■ V7c (100 //m) / 



The coefficients of correlations of the former and latter cases are 0.7 and 0.8, respectively. 
There are 653 and 495 data points for the regions of 7(60 nm)/Ic (100 /xm) < 0.28 and 
> 0.28, respectively. 

We calculate Iq (140 fim) only from 7(60 nm) and Ic (100 /im) for four cases listed 
in Table 1 in addition to the SFD98 method. At first, we use Equation (6) for all data 
points, and call this the single case. As shown in Figure 1, this relation is similar to the 
sub-correlation, and approximately represents all the data points. 

However, the scatter of the data points around the lines of Equations (5) and (6) does 
not reflect the noise included in the DIRBE's individual data points but actual variation. 
The noise level of DIRBE140 is 2 MJy/sr {COBE/BIRBE Explanatory Supplement 1995). 
Actually, the rms noise of DIRBE140 is 2 MJy/sr in the high latitude sky {b > 70°). As 
the lowest intensity is 30 MJy/sr in the present analysis of the Cygnus region, the noise 
does not contribute to the scatter in Figure 1. The noise levels of DIRBE60 and DIRBEIOO 
are negligible compared to the intensities at those bands. Therefore the scatter should be 
actual variation, the fitting with a single line as above introduces an additional error into 
the estimation of temperature. Therefore, we adopted the following method for the latter 
three cases. Each data point in Figure 1 was plotted from the DIRBE data with poor spatial 
resolution. SFD98 determined the dust temperature only from DIRBEIOO and DIRBE140, 
and thus, small-scale variation of the dust temperature was smeared out. On the other 
hand, we can use IRAS60 and IRASlOO with better spatial resolution, and we are aware 
of the reasonably good correlations. Therefore, we initially plotted a single line with the 
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same inclination as that of Equation (5) or (6) just at the data point plotted based on the 
DIRBE. Next, we apply the ratio of IRAS60 and IRAS 100 to the vertical axis of Figure 1 to 
determine the ratio of IRASlOO intensity and the 140 /im intensity. Equation (7) indicates 
this method. 



Ic{UQiim) _ //c(140 /xm)\ // J (60 /xm) \ ^/ / (60 /xm) \ \ 

7c(100 ^xm)IRAs V^c(100 ^m) J ^jj^^^ Vv^c(100 ^J.m)) j^j^g ' V^c(100 nm) ) ^jj^^^J 

(7) 

The power a is -0.75 and -0.44. 

According to the correlation recognized in Figure 1, we adopted three cases as in Table 
1: inclination of Equation (5), inclination of Equation (6), and a combination of both incli- 
nations. We call the case of a = —0.75 the slow case and the case of a = —0.44 the steep 
case. The last case is called the best-fit case that changes the power according to Equation 
(8) for better fitting; 



' ^ -0.44 ((4S/Sl)_ > 0.28) . 

This new method, corresponding to the latter three cases, is based on the fact that the 
scatter around a single fitted line may represent some actual variation and there may still 
remain a correlation with a similar inclination around each data point. It can be interpreted 
that SFD98 adopted a vertical line at each data point for the derivation of the temperature 
in Figure 1, whereas the latter three cases of the present study adopted a tilted line at each 
data point representing the actual correlation found by Hibi et al. (2006). The scatter from 
Equation (5) or (6) produces no error. The DIRBE data were re-gridded from 0.32° pixels 
to 5' pixels, the same as those of IRAS with the sphere-approximation interpolation. As a 
result of this transformation, /c(140 jim) of each pixel can be calculated using Equation (7). 

Figure 2 shows the difference in the methods of calculating the dust temperature among 
five cases: SFD98, single, slow, steep, and best-fit. The best-fit is represented by the combi- 
nation of the steep case and the slow case at /(60 /xm)//c(100 /xm) = 0.28. The difference 
between the present method and SFD98 appears in the difference between the inclinations 
of lines. The deviation of 7(60 //m)//c'(100 jim) — 0.28 of the IRAS data from that of the 
DIRBE data is reflected for Ic (140 /xm)//c(100 jim) = 0.28 that corresponds to the dust 
temperature. Therefore, small-scale variation of the dust temperature should arise between 
the present method and SFD98. 
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In the above mentioned method, the 140 /jtm intensity was calculated from the 60 and 
100 fim intensities of the IRAS data with DIRBE data for the four cases. The Improved 
Reprocessing of the IRAS Survey (IRIS) by Miville & Lagache (2005) was used as the IRAS 
image. The IRIS maps are IPD-subtracted ones, the calibration is rescaled with the DIRBE 
absolute calibration, and the strip patterns are eliminated with a 2-D Fourier filter. 

The IRAS data must be color-corrected as the DIRBE data. The color-correction factors 
were quoted from the IRAS Explanatory Supplement (Beichman et al. 1988). The dust 
temperature was calculated from the observed 100 fim intensity and the derived 140 /xm 
intensity by Equations (6) and (7). The optical depth at 100 iim is written as 

/c(lOO^m) 

rioo^m 5^(100 /xm,T,)- ^ ^ 

Ay is calculated by 

Av = 865.85 X ^loo /.m = 891 x noo ^.m- (10) 

The value of the conversion factor from tioo nm to Ay is quoted from Mathis (1990) by 
assuming Ry = 3.1. This value decreases by 15 % in the case of Ry = 5.0. 

The notation for Ay is summarized in Table 1. 



Table 1: Calculation methods for SFD98 and four new cases. 



Notation Method 



Ay 


(SFD98) 


SFD98 


Ay 


(single) 


A single best-fit hne. Equation (6) 


Ay 


(slow) 


A line with the inclination of Equation (5) at each DIRBE data point 


Ay 


(steep) 


Same as Ay (slow) but with the inclination of Equation (6) 


Ay 


(best) 


Combination of Ay (slow) and Ay (steep) of Equation (8) 



3. Result 

Ay values in the Cygnus region are calculated for the four cases. Figure 3 shows the 
histogram of Ay (best) — Ay (SFD98). The standard deviation is 3.4. Figure 4 shows the 
histograms of Ay (best) — Ay (slow) and Ay (best) — Ay (steep). The standard deviations 
are 1.2 and 0.83, respectively. Figure 5 shows the histogram of Ay (steep) — Ay (single). 
The standard deviation is 1.9. 
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Figure 6 compares the dust temperature maps made by the best-fit case of our study 
with those by SFD98. The dust temperature distribution of the present study shows small- 
scale spatial variation that is not apparent in the map by SFD98. The maximum and the 
standard deviation of the temperature difference between the two maps are 5 K and 0.5 K, 
respectively 

Figure 7 shows the Ay distribution of this area. The fluctuation of both the maps seems 
to have similar angular scales. However, the Ay values of both the maps differ from each 
other. 

To analyze the reason for these differences, the differences in temperature and Ay are 
shown in Figure 8. The sign of the left panel of Figure 8 is reversed compared to that of the 
right panel so as to visually check the dependence because Ay has a negative dependence 
on temperature. It can be seen that the two maps resemble each other remarkably, which 
means that the difference in Ay originates from the difference in the derived temperature. 

4. Discussion 

4.1. Comparison with SFD98 

The methods using the best-fit and steep cases are considered to be more precise than 
the SFD98 method because of the following reason. As described in Section 3, the Ay map 
derived by the best-fit differs from that by SFD98, and the difference can be ascribed to 
the fact that the present temperature map has higher spatial resolution compared with that 
of SFD98. Figure 9 shows the result of comparison between the Ay values of the best-fit 
case and SFD98. The difference, {Ay (best) - Ay (SFD98))/Ay (best), scatters by 21 % 
in 1 sigma. SFD98 employed the IRAS 100 fim intensity to calculate the extinction. They 
removed point sources from the IRAS 100 fim map to smoothen the map with a FWHM = 
3.2' Gaussian profile. The difference in the 100 /xm intensity between our map and that of 
SFD98 is 9 % (1 sigma) in the Cygnus region. This difference can also account for the Ay 
difference between the present study and SFD98 (see Equations (9) and (10) ) in addition to 
the difference between spatial resolution of dust temperature; the scatter seen in Figure 9 is 
affected by this difference. The dust temperature difference in the Ay difference is estimated 
as 19 % by Equation (11); 



AAy{Td) = ^/AAy '^{total) - AAy ^(/(lOO //m). (11) 
AAy{Tci), AAy {total), and AAy {I {100 /im)) is the dust temperature difference at the 
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Av difference, total Ay difference (21 %), and the 100 fxm intensity difference at Ay difference 
(9 %), respectively. On the other hand, the noise of the IRIS data is 0.03 MJy/sr for 
A = 60 fj,m and 0.06 MJy/sr for A = 100 fim (Miville & Lagache 2005). These values 
correspond to less than 1 % for Ay and are negligibly small compared to the difference 
between Ay (best) and Ay (SFD98). 

Figure 10 shows the comparison between Ay (best) and Ay (slow). It can be seen 
that the difference is large compared to the noise of the DIRBE data points, reflecting the 
fact that the power of the slow case does not represent the entire correlation in the Cygnus 
region. Figure 11 shows the comparison between Ay (best) and Ay (steep). The difference 
is 5 % (1 sigma of (Ay (best) — Ay (steep)) /Ay (best)). The difference between Ay (steep) 
and Ay (SFD98) is 20 % (1 sigma of {Ay (steep) - Ay (SFD98))/Ay (steep)). Thus, the 
dust temperature difference in the Ay difference is estimated as 18 % by Equation (11). 
Figure 12 shows the comparison of Ay (steep) and Ay (single). The difference, {Ay (steep) 
— Ay (single))/Ay (steep)), scatters by 21 % in 1 sigma. It is consistent with the hypothesis 
described in Section 2 that the steep, slow, and best-fit cases are more accurate compared 
with the single case. 

4.2. Comparison with the Star Counting Method 

Dobashi et al. (2005) published an Ay map within | 6 | < 40° by the star counting 
method using the DSS images. The spatial resolution is &, similar to the present study. 
Figure 13 shows the result of the comparison of both Ay data in the Cygnus region. The 
Ay (DSS) data is distributed lower than Ay (best) and saturated at Ay (DSS) > 5 mag. 
This indicates that the star counting method with DSS has a limit of around Ay (DSS) ~ 
5 mag. It is apparent that their method is not useful for Ay > 10 mag because the optical 
radiation does not reach us because of heavy extinction by the dust. Dobashi (2009) apphed 
their star counting method to the Two-Micron All Sky Survey (2MASS) to reach Ay ~ 30 
mag. 

4.3. Selecting Valid Correlation 

As shown in Section 4.1, Ay by the best- fit and steep cases are more precise than that 
by SFD98 in the Cygnus region. But these methods would not be effective in high Galactic 
latitudes because Ic (140 pLm) has a lower S/N ratio in those latitudes. Therefore, if we 
choose only one case out of the best-fit, slow, steep, the main-correlation, and the sub- 
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correlation for the entire sky, the sub-correlation is best suited because of its similarity to 
the steep case. 

In the entire sky, there are regions where both the main-correlation and the sub- 
correlation can be seen (e.g., the Galactic plane). There are two methods for determining 
whether the main- or sub-correlation is applicable if the best-fit is not used. The first method 

is the same as that for the Cygnus region. We can make a color-color diagram in any area of 
the sky and choose the better one. The second method involves referring to the radio thermal 
continuum map that represents the bulk of the ionized gas excited by the OB stars. The 
direction with strong radio continuum should be associated with the region where the ISRF 
varies on a small scale because of the strong local heating sources. Therefore, the FIR colors 
in the direction with strong radio continuum can be thought to follow the sub-correlation as 
demonstrated by Hibi et al. (2006). Figure 1 shows that the sub-correlation is more major 
than the main-correlation in the Cygnus region. 

We estimate the systematic difference arising from the selection of either the main- 
correlation or the sub-correlation, since in general we are unaware a priori which of the 
two correlations is applicable to a certain region. We compare Ay calculated by the main- 
correlation and the sub-correlation in Figure 14. The range of the Ay ratio is between 1 and 
2.5 where the range of 7(60 iim) / Ic {100 jim) is between 0.2 and 0.6. 

4.4. Validity of Applying Correlations to the IRAS Map 

Hibi et al. (2006) identified the two groups in the color-color plot of the DIRBE data 
with 42' resolution. Each group has good correlation expressed by Equations (3) and (4). 
These two expressions can be applied to the intensity distribution with a higher spatial 
resolution such as that of the IRAS map. As indicated by Hibi et al. (2006) and Hirashita 
et al. (2007), the main-correlation represents the lines of sight along which the ISRF is 
constant, and the sub-correlation represents the lines of sight along which the ISRF varies 
greatly. Obviously, the lines of sight showing the sub-correlation may also represent the 
direction along which the ISRF varies within the beam. The direction following the main- 
correlation with the DIRBE beam does not include regions having strong radiation (Case 
1). Otherwise, it must follow the sub-correlation. The ISRF does not vary extensively along 
the hue of sight and over the spatial extent of the DIRBE beam. Therefore, all hues of sight 
of the IRAS beam in a DIRBE beam should follow the main-correlation. 

On the other hand, there should be strong ISRF regions in the hue of sight indicating 
the sub-correlation with the DIRBE beam (Case 2). The strong ISRF region is usually 



- 11 - 



compact, so the spatial extent could be smaller than the DIRBE beam. A DIRBE beam can 
be divided into many /i?^S' beams, some indicating the main-correlation (Case 2-1) and the 
others the sub-correlation (Case 2-2). The former does not include the strong ISRF regions, 
and the latter includes the strong radiation field region in the line of sight. 

In Case 2-1, the difference in the derived Ay between the main-correlation and the 
sub-correlation must be small because there is only a small difference between the main- 
correlation and the sub-correlation under weak ISRF. If the sub-correlation is applied, the 
maximum error of Ay is 30 %, where /(60 iim) / Ic{lOO iim) is lower than 0.3 (Fig. 14). 

Finally, we consider Case 2-2. The IRAS data points that include a strong radiation 
field in the DIRBE beam also follow the sub-correlation, because such a strong radiation 
field does not entirely occupy the line of sight. Therefore, there is no error in using the 
sub-correlation in Case 2-2. 

Consequently, even if a higher spatial resolution map (e.g., IRAS) is employed, the Ay 
error does not increase to a large extent. 

4.5. Future Prospect 

An all-sky extinction map can be produced using the present method. In higher Galactic 
latitudes, extinction is generally small and so is the temperature fluctuation. Then the Ay 

difference between the present study and SFD98 would be smaller in the Cygnus region. The 
subtraction of the IPD thermal emission is an important issue for high Galactic latitudes 
or low ecliptic latitudes. The AKARI all-sky diffuse map will be released at 60, 100, 140, 
and 160 ^m. with better spatial resolution in the future. Then, the present study should be 
compared with the AKARI m.&Y>- 

5. Summary 

We have derived the Ay map with 5' resolution in the Cygnus region using the color- 
color correlation between 60, 100, and 140 /im discovered by Hibi et al. (2006). The Ay 
difference between the present study and SFD98 is significantly larger than the uncertainty 
associated with our derivation of the intensity at 140/xm with 5' resolution. The difference 
may occur because the spatial resolution of the present temperature map is ten times higher 
than that of SFD98. Therefore, the present method is more accurate than that of SFD98. 
An all-sky extinction map can be produced by the present method. 
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/c(140 |im)//c(100 |im) 



Fig. 1. — Color-color diagram between 60, 100, and 140 //m in the Cygnus region. The 
horizontal axis is the 140 -100 /im color, while the vertical axis is the 60 -100 //m color. The 
data points arc from DIRBE. The main-correlation, the sub-correlation. Equation (5), and 
Equation (6) are represented by the doted line, dashed line, dashed-dotted line, and solid 
line, respectively. 
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Fig. 2. — Explanation of Ay calculation on the color-color diagram between 60, 100, and 140 
lj,m in the Cygnus region. The horizontal dashed line shows /(60 fim)/Ic (100 fim) = 0.28. 
The circle and square are points arbitrarily selected in regions of /(60 fj,m) /Ic {100 fJ^m) > 
0.28 and < 0.28, respectively. 
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Fig. 3 — Histogram of the Ay difference of Av{best) - Av{SFD98). Bin size is 0.05. 
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Fig. 4. — Histogram of Ay differences. Bin size is 0.05. Filled circles represent the same 
difference as in Fig. 2. 
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Fig. 5. — Histogram of Ay difference of Av{steep) — Av{single). Bin size is 0.05. 
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Fig. 6. — Dust temperature distribution maps in the Cygnus region: (a) present study and 
(b) that by SFD98. The black area in Panel (a) indicates the area that IRAS did not survey. 
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Fig. 7. — Ay distribution maps: (a) present study and (b) that by SFD98. 
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Fig. 8. — Difference maps in temperature (a) and in Ay (b). The sign is reversed for the left 
panel. 
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Fig. 9. — Comparison between the Ay values of the best-fit case and SFD98. 
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Fig. 10. — Comparison between the Ay values of the best-fit case and the slow case. 
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Fig. 11. — Comparison between the Ay values of the best-fit case and the steep case. 
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Fig. 12. — Comparison between the Ay values of the steep case and the single case. 
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Fig. 13. — Comparison between the Ay values of the best-fit case and Dobashi et aL (2005). 
The sohd hne indicates a slope of unity. 
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Fig. 14. — Uncertainty of Ay when the sub-correlation is adopted. The sohd curve shows 
the ratio of Ay by the main-correlation over the sub-correlation. 



